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ABSTRACT 

We present the results of searches for rare B meson decays. The mea- 
surements use all or part of a data sample of about 88 million T(45') 
BB decays collected between 1999 and 2002 with the BABAR detector at 
the PEP-II asymmetric energy B Factory at the Stanford Linear Acceler- 
ator Center. We study a variety of decays dominated by electromagnetic, 
electroweak and gluonic penguin transitions, and report measurements of 
branching fractions and other quantities of interest. 
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1 Introduction 



Measurement j^ ' ^ ' '^ ' ^ ' ^ ' ^ ' ^^ ' ^ ' " ' "' ' ^^ ' ^^ ' ^^ ' ^^ ' ^^ J^ of rare B meson branching fractions have 
been performed using the BABAI^ detector. B decays in which CKM favored am- 
plitudes are suppressed or forbidden are sensitive to penguin amplitudes and hence to 
possible non-Standard Model effects arising from new particles participating in internal 
loops. In addition to probes for new physics, many of these modes are also crucial to the 
full constraint of the "Unitarity Triangle". As the definition implies, rare decays typi- 
cally have branching fractions of less than 10^"^. The present data sample of roughly 
88 million BJ3 pairs allows for measurements or stringent limits on many such modes. 

1.1 Flavor and the Quark Sector of the Standard Model 

The complex CKM^ matrix describes the coupling of the charged weak transition q — > 
W*^q', which is proportional to V^*,. The non-diagonality of this matrix expresses 
the fact that the Weak isospin doublet members (6', s', d!) are states of mixed flavor. 
We can thus view the CKM matrix as the transformation between the mass and flavor 
eigenstates of the quarks 
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The unitarity condition implies that there are four free parameters in this matrix, one 
of which is a phase. It is through this phase that the Standard Model can accommodate 
CP violation. In particular, the orthogonality requirement between the first and third 
columns requires 

KdKT. + v^dV:, + VtdV:, = 0, (2) 

which can be expressed geometrically as the so called "Unitarity Triangle" shown in 
FigurelU Information about each side and angle is accessible through a variety of mea- 
surements in the B meson system. The angles are measured through time-dependent 
decay rate asymmetries, and the sides via direct or indirect measurements of the CKM 
matrix elements. Measurement of all the components of the Unitarity Triangle over- 
constrains the triangle, and thus provides a test of the Standard Model (SM). 




Fig. 1. The Unitarity Triangle determined from the orthogonality of the first and third 
columns of the CKM matrix. Also shown are B meson processes which yield informa- 
tion about each side and angle. 

2 The BaBar Detector 

A detailed description of the BABAR detector can be found elsewhere.'^ Charged parti- 
cle momenta are measured in a tracking system that consists of a 5-layer double-sided 
silicon micro-strip vertex tracker (SVT) and a 40-layer drift chamber (DCH) filled with 
an (80:20) mixture of helium and isobutane. The tracking volume is contained within 
the 1.5T magnetic field of a superconducting solenoid. The combined track momentum 
resolution is o-p^/pT = 0.13% x px + 0.45%. The primary charged hadron identifia- 
tion device is a detector of internally reflected Cerenkov radiation (DIRC). The typical 
separation of kaons and pions due to their measured Cerenkov angle 9c varies from 8cr 
at 2 GeV/c to 2.5a at 4 GeV/c, where a is the average 9c resolution. Specific ioniza- 
tion energy loss (dE/dx) measurements in the DCH and SVT also contribute to charged 
hadron identification for particle momenta less than 0.7 GeV/c. Photons are detected in 
an electromagnetic calorimeter (EMC) consisting of 6580 Thallium doped Csl crystals 
arranged in barrel and forward end-cap sub-detectors. The 7r° mass resolution in on 
average about 7 MeV/c^. Muons and long-lived neutral hadrons are detected within the 
instrumentation of the solenoid flux return (IFR) which consists of alternating layers of 
iron and resistive plate chambers. 



3 Common Analysis Features 



3.1 Data sample 

The analyses described here use all or part of a data sample consisting of approximately 
88 million pairs of T(4S') BB decays, corresponding to a detector exposure of about 
81 fb^^ An additional sample of 9.6 fb"^ taken about 40 MeV below the peak of the 
T(4S') resonance ("off-resonance") is used by many analyses to study e+e^ — qq 
"continuum" backgrounds. 

3.2 B Meson Reconstruction 

B mesons produced from T{AS) decays are identified via their unique kinematics. Be- 
cause the mass of the B meson pair is nearly that of the T(45'), they are produced 
nearly at rest (p^ ^ 325 MeV/c). Use of the beam energy in constraining the kinemat- 
ics serves to reduce the resolution of theses variables. 
The conservation of energy can be expressed as: 

AE = E*^- (3) 

where £'beam is the single beam energy in the center-of-mass (CM) frame. E]^ is the 
measured energy of the B candidate in the CM. Correctly reconstructed B candidates 
have AE distributed around zero with a resolution ranging from 15 to 80 MeV. The 
energy resolution of the B decay products dominates the resolution of this variable. 
Continuum background in this variable is well described by a monotonically decreasing 
low order polynomial. Figure |2l shows the AE distribution for a typical rare mode after 
all selection criteria have been applied (except that on AE). 
We express momentum conservation as: 

rnES = ^Je^I^-p"^. (4) 

Here niEs is the "beam-energy substituted mass", with the S-candidate momentum 
in the CM. Correctly reconstructed B candidates have niEs equal to the B meson mass, 
with a resolution of about 2.5-3.0 MeV/c^, which is dominated by the beam energy 
spread. The continuum background shape in tres is parameterized by a threshold 
functior^ with a fixed endpoint given by the average beam energy. Figure 13 shows the 
Mes distribution for a typical rare mode after all other selections have been applied. 
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Fig. 2. AE distribution for charged B decays to three charged kaons. All candidate se- 
lection criteria have been applied except that on AE. The solid line shows the expected 
continuum background level. 

In addition to the kinematics of the B meson, signal events are selected by mak- 
ing requirements on the decay products. B daughter resonances are required to have 
invariant masses within a restricted range typically determined by resolution and the 
need to leave sufficient sideband to determine background levels. Particle identifica- 
tion requirements are made to select some particles and veto sources of background. 



3.3 Background Suppression 

All rare analyses suffer from substantial backgrounds, and a variety of techniques are 
employed to reduce this to manageable levels. In general, backgrounds from other B 
decays are small. Decays resulting from CKM favored b c transitions have heavier 
daughters and higher multiplicity final states than do CKM suppressed decays. In order 
to wrongly reconstruct such a decay as a rare signal, one must typically lose a particle 
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Fig. 3. rrtEs distribution for charged B decays to three charged kaons. All candidate 
selection criteria have been applied except that onmEs- The solid line shows the ex- 
pected continuum background level. 

from the true B decay, resulting is a substantial shift in the candidate's l\E. The only 
exception to this is for rare modes with high final state multiplicities, in which indi- 
cations of some background have been observed. While B decays from other CKM 
suppressed transitions have similar kinematics and multiplicities to that of a rare signal, 
such modes are rare themselves, and require only limited suppression in most analyses. 
Where B backgrounds are present, they typically populate the sidebands of the AE' 
distribution, but have tails that reach into the signal region as illustrated in Figure HI 

For all the modes discussed here, the primary background is due to random particle 
combinations arising from continuum quark-antiquark production. Although the prob- 
ability for any given continuum event to satisfy a signal selection is quite small, the 
numbers favor the continuum. The total production cross section for light quarks (in- 
cluding charm) under the T(4S') is about 3.5 nb, but only 1 nb for the T(4S') itself. For 
a mode with an expected branching fraction of order 10^ this means that continuum 
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Fig. 4. Typical niEs vs distribution after event selection. The signal populates the 
region around = 5.280 and = 0. Continuum background populates the entire 
plane. B background populates the A£J sideband. 

events are produced at a rate well in excess of 10^ times that of the signal. 

In order to control continuum backgrounds, one typically exploits the fact that while 
the B meson pairs are produced near threshold in T(45') decays, the light and charm 
quark pairs which comprise the continuum are produced with a great deal of excess 
energy. The result is that for true B meson decays, the decay products entering the de- 
tector are distributed isotropically in the CM, while the continuum background exhibits 
a "jet-like" topology, with a strong correlation between the B candidate decay and jet 
axes. 

The first topological variable typically employed is the angle Or between the thrust 
axes of the B candidate and the remaining particles in the event. The sphericity axes 
may be used almost interchangeably. The distribution of \cos9t\ is nearly uniform for 
true B mesons, but is strongly peaked near 1 for continuum background as is illustrated 
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Fig. 5. Signal and background distributions of \cos6Thrust\ for a typical B meson decay. 
The signal distribution (solid histogram) is uniform, reflection the random orientation 
between the B candidate thrust axis and the thrust axis of the rest of the event. The 
background (open points) is strongly peaked at one due to the strong correlation be- 
tween these two axes in fake B candidates arising from random particle combinations 
in the continuum. 

in Figure |5] If additional background rejection is required, one may consider the re- 
maining event shape information, such as the angles between the B thrust and decay 
axes and the beam as well as the angular energy flow in the event, and combine it into 
an optimized quantity using a neutral network or a Fisher discriminant. An example of 
the the separation power of a Fisher discriminant after a thrust cut has been made can 
be found in Figure |6l 

3.4 Signal Extraction 

All rare analyses at BABAR are performed "blind", meaning that signal yields are hidden 
from the analyzer until the analysis has been peer reviewed and determined to be in a 
final form. These steps are taken to avoid experimenter's bias. 

There are two primary methods in which signal event yields are extracted; the event 
counting analysis, and the maximum likelihood fit. Detection efficiencies determined 
from signal Monte Carlo simulations and data control samples are used to convert yields 
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Fig. 6. Distribution of Fisher discriminant for data control mode (solid points), control 
mode signal Monte Carlo (solid histogram), continuum data (open points) and contin- 
uum Monte Carlo (dashed histogram) after a cut on \cos6Thrust\- The Fisher and Thrust 
angle are strongly correlated, thus the separation will depend strongly on the thrust cut 
made. 

into branching fraction measurements. Equal production of charged and neutral B's 
from the T (45") decay are assumed throughout. 

In the event counting analysis, a set of selection criteria are defined to select a signal 
region of the parameter space. The criteria are optimized with respect to expected signal 
and background yields to produce a measurement of the greatest possible statistical 
significance. Systematic uncertainties may figure into this process, but they are usually 
negligible for rare modes. The selection are applied to the data, and the population of 
signal region is counted. An estimated background is subtracted to determine the signal 
event yield. The background yield is typically determined by measuring the density of 
events in a sideband region and projecting that density into the signal region. 

Li the maximum likelihood fit method signal yields are determined by an unbinned 
extended maximum likelihood fit to a set of observables. These typically include uies, 
AE, event shape variables and, where appropriate, B daughter resonance invariant 
masses and particle identification information. The probability Vi{xj] Si) for a given 
hypothesis i is the product of probability density functions (PDFs) for each of the vari- 
ables Xj = {mEs, AE, Fisher, ...) given the set of parameters ai. The hypotheses i are 



signal, continuum background and sometimes B background for each final state in the 
fit. The likelihood function is given by a product over all events and the signal and 
background components: 

The Hi are the numbers of events for each hypothesis. The values of the yields (and any 
free parameters in the PDFs) are taken as those which maximize the likelihood function. 
Unit change in —2lnC defines the one standard deviation statistical uncertainties on the 
free parameters in the fit. The statistical significance of the signal yield is determined 
from the change in —2lnC when the signal yield is forced to zero. If no statistically 
significant signal is found (more than 4 standard deviations), a 90% confidence level 
upper limit may be obtained by requiring: 

C{n)dn 

Full and toy Monte Carlo simulations are used to verify that the fit is unbiased. 

The accuracy with which the PDFs describe the data is of utmost importance in the 
likelihood fit. Background PDFs are determined by fits to off -resonance and sideband 
data. Signal PDFs are determined primarily from signal Monte Carlo simulations, but 
ultimately rely on data control samples to verify their validity. 



4 Electromagnetic Penguins 

Electromagnetic penguins consist of the class of amplitudes in which an external photon 
is emitted by one of the virtual particles participating in the loop through which the b 
s{d) transition proceeds. This is illustrated in Figure IT] Such diagrams are relatively 
clean from a theoretical perspective, and a variety of information can be gather from 
measurements of decays dominated by these amplitudes. The decay B K*'y was the 
first penguin to be observed.!^ Measurements of its branching fraction provides a test 
of QCD, and direct CP violation in this mode would be an indication of new physics. 
The decay rate ratio of B ^ p'j to B ^ K*'y is sensitive to the ratio of \^\- The 
photon energy spectrum from measurements of 6 — > 57 provides information on the 
mass and Fermi motion of the b quark within the B meson. 

In the analysis of these modes, in each case there is a requirement of a high energy 
isolated photon. The calorimeter cluster is required to have a profile consistent with 




Fig. 7. Feynman diagram for an electromagnetic penguin amplitude. 

an electromagnetic shower, and the candidate photon must not be consistent with hav- 
ing originated from a vr" or 77 decay. Further details and results of each analysis are 
presented below. 

4.1 Measurement of B ^ K*j 

The analysis of B ^ K*'y^ has been performed on a data sample corresponding to 
approximately 22 million BB pairs recorded in 1999-2000. The K*'j final state is 
reconstructed in all four K* decay modes. Stringent identification requirements are 
placed on charged kaons. Invariant mass requirements are placed on both K* and Kg 
candidates. The Kg is also required to have a decay vertex displaced from the e~^e~ in- 
teraction point. Since the B meson is a pseudoscalar, angular momentum conservation 
requires that the K* is polarized. The absolute value of the cosine of the K* helicity 
angle is required to be less than 0.75. Continuum background is suppressed with cuts 
on the absolute values of the cosines of the thrust and B flight angles, both of 0.80. 

After cutting on AE, the signal yield is determined from an unbinned maximum 
likelihood fit to the ities distribution, shown for each K* decay mode in Figure |H1 
Branching fraction and direct CP asymmetry results are shown in Table [Tpi' 





B{B^ K*^^) 


BiB+ K*+-i) 


Acp 


Theor>22 23 24^ 


7.5 ±3.0 


7.5 ±3.0 


Acp < 0.005 


BABAR 


4.23 ±0.40 ±0.22 


3.83 ±0.62 ±0.22 


-0.17 < Acp < 0.08 
@90% CL 



Table 1. Results of the branching fraction and direct asymmetry analysis of i? — > 
K*'y. In each result, the first uncertainty is statistical, the second systematic. Branching 
fractions are in units of 10~^. 
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Fig. 8. The niEs distribution for B — > K*7 in each of the K* decay modes. All selec- 
tion criteria have been applied. The solid curve is the combined signal and background 
PDF shape, the dashed is the background only. 



4.2 Search for B ^ and B ^ 

The analysis of the p7 and cwj final state^^is significantly more challenging than that of 
K*^. The predicted branching fractions are about 50 times smaller than for K*'j. Both 
the p and the to have significantly more background under their peaks than does the K*, 
and the rho is much broader. In addition to continuum background, these modes also 
potentially suffer cross-feed background from K*'j, other b ^ s'j processes and from 
B p7r°. 

A neutral network containing information from event shape. At and flavor tagging 
is used to control continuum background. K*7 feed-across is vetoed using particle 
identification. After these selection criteria are applied, the signal yield for each final 
state is extracted using a unbinned maximum likelihood fit to mEs, ^E, and the p/cu 



invariant mass. Studies of generic BB Monte Carlo show that the expected B back- 
ground is quite small, so the fit includes components only for signal and continuum 
background. Background from B decays is considered as a systematic uncertainty. 

The results of this analysis applied to a sample of 84 million BB pairs can be found 
in table|2| If isospin symmetry is assumed, all three modes can be combined to produce 
an upper limit of B{B p'j) < 1.9 x 10"'' @ 90% CL. This resuk can be used to 
place on upper limit on CKM parameters |^| < 0.036 @ 90% CL. A discussion of 
theoretical errors can be found in Ali and Parkhomonko.*^ 





B{B^ ^ p07) 


B{B+ ^ p+7) 


B{B^ CU7) 


Theor>^4 


0.5-0.75 


0.8 - 1.5 


0.5-0.75 


BABAR 


< 1.4 


< 2.3 


< 1.2 



Table 2. Results of the branching fraction analysis of B — p7 and uj^y. Branching 
ratios are in units of 10"^. Upper limits are at 90% CL. 
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4.3 Semi-Inclusive Measurement of 6 ^ 57 

This ana lysiPof 22 million 55 pairs is a study of a collection of exclusive final states 
with a kaon plus up to four pions, no more than one of which may be neutral. Because 
— 57 is a two-body decay process, the photon energy in the B rest frame is related 
to the recoil hadronic mass, Muad' 

_Ml-Ml 

Fits to the measured spectra of both of these quantities can be used to determine the total 
branching ratio for B — > Xg'y^ In addition to constraining new physics contributions 
to the underlying amplitude, parameters associated with heavy quark effective theory 
(HQET) are also extracted in the analysis. These parameters are critical to reducing 
theory errors in the extraction of Vub and Vcb- 

Measured branching fractions as a function of Muad and can be found in Fig- 
ure |9l Analysis of these spectra yield results: 

A = 0.37 ± 0.09(stat) ± 0.07(s?/st) ± 0.10(moc/e/) GeV/c^ 

mb = 4.79 ± Om{stat) ± 0.10{syst) ± 0.10{model) GeV/c^ 

Ai = -0.24l°:°3 (stat) ± 0.02{syst) tofi{model) GeV/c^ 

B{b 57) = 4.3 ± 0.5{stat) ± 0.8{syst) ± 1.3{model) x 10"^ (8) 
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Fig. 9. Hadronic recoil mass and photon energy spectra for the semi-inclusive analysis 
of 6 ^ 57. Errors are statistical only. 

4.4 Fully-Inclusive Measurement of 6 ^ 57 

Much of the uncertainty in the semi-inclusive 6 — > 57 analysis arises from theoretical 
errors. HQET implies a duality between the quark and hadron level of an interaction, 
which implies that parton level rate for 6 ^ 57 is the same as the inclusive rate for 
B Xg'j. These two issues motivate the fully inclusive analysis technique. 

This analysiPis performed on a sample of 60 million BB pairs. Photons in the 
range 1.5 < E* < 3.5 GeV are analyzed. These photons are required to meet the selec- 
tion criteria described above. To suppress continuum background the event is required 
to have a lepton flavor tag which strongly selects true BB decays. Such a selection 
induces no model dependency into the analysis as it only applies to the "other" B in 
the decay. In addition to particle identification criteria, fake leptons are further rejected 
by requiring large missing energy in the event, which is normally associated with semi- 
leptonic b ^ c transitions. To obtain additional discrimination, the angular separation 
between the lepton and the photon is required not to be small. Event topology in the 
form of the Fox-Wolfram moments in also employed to reduce background from the 
continuum. Backgrounds are estimated using off-resonance data and BB Monte Carlo. 

The spectra for on-resonance data and the predicted background are shown in 
Figure fTOl The photon energy range 2.1 < E* < 2.7 GeV is considered to reduce 
model dependencies. The branching fraction for B Xs'j is measured in this region 
and then extrapolated to the full spectrum: 

B{B X,7) = 3.88 ± 0.36(stat) ± 0. 37 {syst)totl{model) x 10"^ (9) 




E, (GeV) 

Fig. 10. The E* distribution of on-resonance data (solid points) and background expec- 
tations for the fully inclusive analysis of B ^ Xs'j. Errors are statistical only. 



5 Electroweak Penguins 

As the name suggests, Electroweak penguins are amplitudes that proceed via loops 
involving photons, W or Z bosons. Such processes are strongly suppressed in the 
Standard Model, and as such, are excellent windows onto potential new physics. Well 
controlled theoretical uncertainties aid in this sensitivity. We discuss the analysis of 
four such final states here: K '■*k+ i-^i+i-^KuuP and 77.EI 



5.1 Measurement of B ^ K'^*k+i- 

The flavor-changing neutral current decays B Ki^i^ and B K*i'^£^ have 
predicted branching fractions on the order 10^^ — 10"''.^^ The leading diagrams for 
this decay are electroweak penguin and box diagrams, and can be found in Figure [TT] 

The decay rate for B K^*H'^i~ is rather sensitive to the presence of new physics. 
In particular, certain extensions to the SM can vary the rate by more than a factor of two. 
In addition to the decay rate, kinematic distributions accessible with higher statistics, 
such as the boson distribution (mj^) and the forward-backward asymmetry in the K* 
channel are of considerable interest as they are also quite sensitive to non-SM physics, 
and are less model dependent than the overall rate. 
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Fig. 11. Leading Electroweak penguin and box diagrams for the decay B — > K'^*H'^l~ . 

The experimental challenge in this analysis is to control the various sources of back- 
ground. Background from B charmonium decays which have the same final state 
particles are control by vetoing regions in the /\E vs niu plane. Continuum background 
is reduced using a Fisher discriminant which in addition to event shape information 
includes information on the K^ invariant mass, which serves to veto D Kiv. Com- 
binatorics from semi-leptonic B decays are rejected using a B-likelihood built from the 
missing energy in the event, vertex information, and the B production angle. Finally, 
peaking backgrounds from particle mis-identification are reduced by vetoing the K'^*^tx 
mass in the region of the D mass. 

After background rejection and particle identification criteria are applied, the signal 
is extracted with a likelihood fit to rriEs and AE. The results of this analysis on a 
sample of 88.4M BB pairs are: 

B{B K£+r) = (0.78l°J^(stat)+°;}^(s?/st)) x 10"^ (10) 

with a significance (including systematics) of 4.4cr, and 

B{B K*e+e-) = {l.GStofsistat) ± 0.28{syst)) x 10"^ (11) 

with a significance of 2.8cr. Since the K* result is not significant, we report a 90% CL 
upper limit: 

B{B K*tr) < 3.0 X 10"^ (12) 
Combined projections of niEs and are shown in Figure fT^ 
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Fig. 12. Projections of niEs and for the combined K^*H^£ final states. The solid 
line is the combined signal plus background PDF, the histogram are the data. 



5.2 Search for B 







£+t 



The decay of a i? meson to a pair of leptons is highly suppressed within the SM by 
factors resulting from CKM, internal quark annihilation and helicity. Leading diagrams 
are shown in Figure[T3l Within the SM, predicted branching fractions are 1.9 x 10^^^ 
and 8.0 x 10^^^ for the e+e~ and /i^/i^ channels respectively .'^^ The efi channel is 
forbidden by lepton number conservation. New physics can significantly alter these 
predictions.!^ 

In this analysis, the primary sources of background are from real lepton production 
from continuum cc decays, pions which are mis-identified as muons, and two photon 
processes. Continuum background is suppressed using the thrust magnitude, and the 
angle between the thrust axes of the B candidate and the rest of the event. A track mul- 
tiplicity cut serves to reject two photon processes. The signal is selected by requiring 
two high momentum leptons of opposite charge and good vertex information. Particle 
identification requirements for both leptons are made. The signal yield is determined by 



Fig. 13. Leading diagrams for the decay B — > i^i 



counting events in a signal region of m es and AE, and subtracting an estimated back- 
ground determined from the scaled population of the niEs vs AE plane. The results of 
this analysis applied to a sample of approximately 60 million BB pairs are presented 
in Table |3l 





Ngsb 


NsigBox 


Nbg 


90% CL Upper Limit 


B{B^ e+e-) 


25 


1 


0.60 ±0.24 


3.3 X 10-^ 


B{B^ fi+fi-) 


26 





0.49 ±0.19 


2.7 X 10-^ 


B{B^ fi+fi-) 


26 





0.49 ±0.19 


2.7 X 10-^ 



Table 3. Results of the search for B £~^£~. Ngsb is the population of the tues 
AE sideband and NsigBox the population of the signal region after all selection criteria 
have been applied. Nbg is the expected background in the signal region based on the 
sideband population. 

5.3 Search for 5+ K+uu 

Within the SM, the decay h suu is a pure electroweak flavor changing neutral cur- 
rent. The final state is nearly free of strong interaction uncertainties, and hence the 
theoretical errors associated with this decay are small. While the inclusive analysis is 
not currently feasible, it is possible to search for the exclusive decay B^ K^vu. 



Summing over all neutrino species, the SM prediction for this branching fraction^^ is 

B{B+ K+vv) = 3.8iJ:^ X 10-^ (13) 

The presence of two neutrinos in the final state makes this analysis difficult, as there 
are no kinematic constraints which may be applied to the signal B. Instead, the strategy 
is to fully reconstruct the other B from the T{AS) decay, and compare the remaining 
particles in the event with the signature expected from the signal. The "tag" B is re- 
quired to be fully reconstructed as either B^ D^t^u or B^ D*^£~u. The D° is 
reconstructed in the K~7r'^, K^n^n^n^ and K^n^n^ modes, which results in a total 
of about 0.5% of all charged 5s being reconstructed as tags. To select signal events, 
a high momentum charged kaon is required in the recoil of the tagged B. Additional 
requirements are made on the neutral energy in the recoil and on the angle between the 
kaon and the tag side lepton. 

10 



S 8 



6 



sa 4 

Q 2 
^ 
-2 
-4 
-6 
-8 
-10 



_l 1 1 1 1 1 1 1 1 

_• 

• 

• 


1 1 1 1 1 1 1 1 
• 


, , 1 , ^1 , , , 1 , , , , 1 , , , 1 1 , , 1 1 , , i_ 

• 

. BABAR 

• • - 

• 

• 


• 

• 


• 

• 


• 

• -• 

' • *. 

• / . •. 

. ••• . ' 

• .- 


• 

• 

data 

^ < . . 1 . . . < 1 . . . . 1 < . . 


• • 
• • 

• 

• • * 
• • • . - 

. . < . 1 < . . . 1 . . . . 1 . . . <j . < . . 1 < . . r 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 2.25 2.5 

Remaining neutral energy (GeV) 



Fig. 14. Distribution of K^uv candidates in the plane defined by the electromagnetic 
energy not from the tag B, Eieft, and the resolution scaled difference between the 
measured and mean measured D mass. The signal region is the bounded area at the left 
of the plot. 



Events are counted in a signal region in the plane defined by the electromagnetic 
energy in the tag B recoil, and the difference between the reconstructed and mean fitted 
D mass, scaled by the fitted D mass resolution. The expected background, determined 
by scaling the sideband population into the signal region, is subtracted from the signal 
region population to determine the signal yield. This is illustrated in Figure [14| In a 
sample of 60 million BB pairs, the expected background in the signal region is 2.2 
events, and there are two events observed. The 90% confidence level upper limit on the 
branching fraction, including systematics, is 

^ K+uu) < 9.4 X 10~l (14) 

5.4 Search for 77 

The decay B^ 77 is an example of electroweak annihilation. The SM expectation 
for this decay is small, with predictions ranging from 0.1 to 2.3 x 10^*^.^° As with the 
other modes discussed in this section, physics beyond the SM can result in significant 
enhancements to this rate.l^ 

In this analysis, selection criteria are placed on the ratio of the 2nd to 0th Fox- 
Wolfram moments, the cosine of the angle between one of the photons (chosen at ran- 
dom) and the thrust axis of the rest of the event, and the B production angle to suppress 
continuum background. Selected photons are required not to be consistent with having 
come from a vr^ or 77 decay. The signal yield is determined by counting events in a 
signal region of the plane defined by tues and AE, and subtracting the expected back- 
ground determined by scaling the the sideband population into the signal region. The 
result for a sample of 22 million BB pairs is 

^ 77) < 1.7 X IQ-^ (15) 

at the 90% confidence level, including systematic uncertainties. 

6 Gluonic Penguins (Charmless Hadronic B Decays) 

Charmless hadronic B decays proceed through a combination of CKM suppressed tree 
(b u) and gluonic penguin (6 d, s) amplitudes. There are about 70 possible com- 
binations of two-body decays in the lowest pseudoscalar and vector nonets. These may 
be further broken into two groups; two-body decays in which both B daughters are 



kaons or pions, and quasi-two-body decays in which at least one of the B daughters is 
a short-lived resonance. The two-body modes can be analyzed for information on the 
CP phases a and 7, and have been found to have significant penguin contributions in 
addition to CKM allowed tree amplitudes. Several of the quasi-two-body modes are 
sensitive to the CP phase p. In addition to yielding information about the Unitarity Tri- 
angle, decays in which penguin amplitudes are dominant are sensitive to new physics. 
Our study of three-body B decays has thus far been limited to combinations of three 
charged kaons or pions. 

All of these modes share some common features. The primary source of back- 
ground is random particle combinations in the continuum, although modes with large 
final state multiplicities or significant neutral energy may suffer from non-negligible 
BB backgrounds. All the final states are ultimately composed of high momentum 
kaons and pions, so the ability to distinguish been these particles at high momenta is 
crucial. 

6.1 Two-Body Decays 

Two body B decays to kaons and pions are sensitive to the angle a of the Unitarity 
Triangle through the time-dependent CP violating asymmetry in the decay B tt+tt^ 
and to the the angle 7 through branching fractions and direct CP- violating asymmetries 
of decays to various tttt and Kn final states. Because there are substantial penguin 
amplitudes which contribute to the tt+tt final state in addition to the tree amplitude, 
the time-dependent asymmetry in that mode does not directly measure a. An isospin 
analysis of the rates for all the B ^ nn decays is required to fully unfold the effects 
of the penguin contributions and determine the relationship between what is measured 
from the tt+tt^ analysis (ag//) and a. Interference between penguin and tree ampli- 
tudes may also lead to substantial direct (time-independent) CP asymmetries in the Kn 
final states. 

In each mode, the signal is extracted using an unbinned extended maximum likeli- 
hood fit, using the tues, ^E, a Fisher discriminant, and where appropriate, Cerenkov 
angle residuals. Groups of related decays are fit simultaneously. For example, the 
TT+TT^, K^n^ and K^K^ yields are determined from a single fit. In these cases, /\E 
and the Cerenkov angle residuals separate the signal modes from each other. Branching 
fraction results for all two-body modes based on a sample of 88 million BB pairs can 
be found in Table |4| 



Decay 


Nsignal 


B X 10-*^ 


Acp 




157 ± 19 


4.7 ± 0.6 ± 0.2 






589 ± 30 


17 9±09±07 


-0 102 ± 050 ± 016 




1 ± 8 


< 0.6 




B+ 7r+7r° 


125 ± 22 


55±10±06 


-0 03 ± 18 ± 02 

KJ B KJ \J — 1 — KJ * -L \^ — 1 — KJ * KJ ^ 


B+ ^ K+7i° 


239 ± 22 


12.8 ± 1.2 ± 1.0 


-0.09 ± 0.09 ± 0.01 


B+ ^ i^o^o 


86± 13 


10.4 ± 1.5 ±0.8 


0.03 ±0.36 ±0.09 




23± 10 


<3.6 (1.6i°il°J) 




B+ ^ K^T,+ 


172 ± 17 


17.5 ± 1.8 ± 1.3 


-0.17 ±0.10 ±0.02 




< 10 


< 1.3 





Table 4. Results of two-body branching fraction analyses.'^'^^'^^^ The 7r°7r° result has 
a statistical significance of 2.5cr. The results for K°7r+ and K^K^ are based on 60 
million BB pairs. Upper limits are at the 90% confidence level. 

Despite no central measurement of the tt^tt^ final state, it is still possible to place 
limits on the relationship between the measured parameter tte// and the Unitarity Trian- 
gle parameter a. Using the bound of Grossman and QuiniP^and our measured values, 
we set an upper limit of |ae// — a| < 51° at 90% CL. 

6.2 Quasi-Two-Body Decays 

Quasi-two-body decays proceed through resonant intermediate states. The analysis 
of such modes is very similar to true two-body decays, but there are additional vari- 
ables that provide separation between the signal and background, such as the reso- 
nance invariant mass and polarization (if the final state is a pseudoscaler-vector com- 
bination). We present the analyses of three groups of related quasi-two-body decays; 
B (j)K^*\^^ B cj/PEland B r]^')K^*)d5, 

6.2.1 B (l)K^*^ 

The decay b — > sss is CKM forbidden, thus the decay B — > (pK^*^ is a nearly pure 
gluonic penguin, as shown in Figure [T51 New physics might not only manifest itself 
as a deviation from the SM prediction for the decay rate, but since the mode (pKg is 
the strange analog to J/iljKg, the CP phase one measures in a time-dependent analysis 



s 
s 

u, d 

Fig. 15. Leading diagram for the decay B — > (j)K^*\ 

could be altered from it's SM value of (3. 

The signal yield in each of these modes is determined from an extended unbinned 
maximum likelihood fit to niEs, ^E, a Fisher discriminant and the (p {K~^K~) invari- 
ant mass. For the (j)Kg, (j)K^ and 07r+ final states, the <p polarization is included in the 
fit, as is the Cerenkov angle residual for the charged states. For the (pK* final states, 
the K* invariant mass is included in the fit. Significant signals are observed for both 
charged and neutral B decays to K and K* final states. The results of this analysis on 
a sample of 60 million BB pairs are 




B{B^ - 


^<PK+) = 


B{B+ - 


-.<PK') = 


B{B+ - 


>^K*+) = 


B{B+ - 


= 


B{B+ 





(9.2 ±1.0 ±0.8) X 10"*' 
(8.7+};^ ±0.9) X 10"^ 

(8.7^1? ±1.1) X 10-6 

0.56 X 10-*^ @ 90% CL. (16) 

07r+, which is both CKM and color 



A stringent limit is also placed on the decay B 
suppressed. 



6.2.2 B ^ uh {h = K, vr) 

B decays involving an omega and either a kaon or a pion proceed through a mixture of 
CKM suppressed b ^ u tree and CKM forbidden b ^ d,s penguin amplitudes. The 
analysis method is identical to that described for (f)K^*\ The results in Table|5]are based 
on a sample of 22 million BB, except for the ujKg analysis, which was performed on 
a sample of 60 million BB pairs, and is a first observation. 



Final State 


N Signal 


Sia) 


B X 10-6 


uK+ 


6.4t|i 


1.3 


<4 (1.4li;|;±0.3) 


ujK^ 




6.6 


5.9tl:^±0.9 


un^ 


27.6t?:? 


4.9 


6.6lfi±0.7 




-0.911° 




< 3 (-0.3 ±1.1 ±0.3) 



Table 5. Results of the branching fraction analysis of i? ^ uh. S is the statistic 
significance of the result. Upper limits are at the 90% confidence level. 

6.2.3 B r^W/sT^*) 

B decays to r] and i]' with a kaon or K* proceed predominantly through penguins, 
although there is some b ^ u tree contribution as well. The decays B rj'K and 
B TjK* were the first gluonic penguins to be observed^ and the rates are much 
larger than initially expected. The best present conjectur^is that the tree and penguin 
amplitudes interfere in such a way as to enhance "q'K and rjK* but suppress rj'K* and 
r]K. Because of it's relatively large rate and nearly pure penguin content, B v'^s 
is also of considerable interest for measurements of time-dependent CP asymmetries, 
which within the SM should probe the angle /3. 

Signals for these modes are extracted as described above for cu and </>. The r]' is 
reconstructed in two decay chains; ?7(77)7r+7r- and p°7. The ?/ is reconstructed as 
^ 77 and 7] n^n^n^. The results of these analyses are displayed in Tabled The 
data samples used for the rj and i]' analyses are 22 and 60 million BB pairs respectively. 



Final State 


-^Signal 


B X 10-6 


7]'K+ 


445 ± 26 


67±5±5 




135 ± 15 


46 ± 6 ± 4 




5.2 ±3.4 


< 13 (4.0t|^ ± 1.0) 


7]K+ 


12.9 ±5.7 


< 6.4 (3.811:^ ±0.2) 




8.0 ±5.9 


< 5.2 (2.2ti| ±0.1) 


r]K° 


5.7 ±3.3 


< 12 {G.Otli ± 0.4) 


r]K*^ 


20.5 ±6.3 


19.8li| ± 1-5 


r]K*+ 


14.3 ±6.6 


22.11^^2^ ± 3.2 



Table 6. Results of the branching fraction analyses of B ^ ri'''^K^*\ Upper limits are 
at the 90% confidence level. 



6.3 Three-Body Decays 

We describe here the analysis of B'^ —* h~ h'^^ where h is either a charged kaon 
or pion. An event counting analysis is performed over the full three particle dalitz 
plot. All final states are measured simultaneously, and unfolded to obtain branching 
fractions for each combination. Continuum background is suppressed using the thrust 
angle and a Fisher discriminant. In addition to continuum background, the open nature 
of the dalitz plot also admits background in some regions from i?+ J/^K^ and 
5+ — i> Dn^ / DK^ . These regions of the dalitz plot are vetoed. Charged particle 
identification is crucial to this analysis, and along with tracking, is the primary source 
of systematic uncertainty. Figure [T^ shows the dalitz plots for 5+ K^K^K^ and 
— i> K^n^n^ . Results of this analysis on a sample of 56 million BB pairs are 

B{B^ ^ -K^Ti^-K^) < 15 X 10"^ @ 90% CL (8.5 ±4.0 ±3.6) 

B{B^ K^TT^TT^) = (59.2 ± A.7{stat) ± A.9{sys)) x 10"^ 

B{B^ ^ K^K^n^) < 7 X 10"^ @ 90% CL (2.1 ±2.9 ±2.0) 

B{B^ K^K^K^) = (34.7 ± 2.0{stat) ± l.8{sys)) x 10"^ (17) 
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Fig. 16. Unbinned dalitz plots for K^n^Ti^ (left) and K^K~K~^ (right) for events in 
the signal region. No efficiency corrections have been applied to the dalitz plots, and 
the charm contributions have not been removed. 



7 Conclusion and Outlook 



We have presented a number of results for rare B meson decays using all or part of a 
sample of approximately 88 million BB pairs collected by the BABAR detector. Updates 
of many of these analyses to the full data set are in progress. These results represent 
only a part of the spectrum of possible measurements of rare decays. The larger data 
sets that will be available in the coming years will allow us to more fully exploit rare 
decays to test the self consistency of the flavor sector of the Standard Model, and will 
perhaps offer the first glimpse of new physics which lies beyond. 
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